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| Abstract 

00 ■ 

f — . We investigate the relation between the Riesz and the Baez-Duarte cri- 

terion for the Riemann Hypothesis. In particular we present the relation 
between the function R(x) appearing in the Riesz criterion and the se- 
quence Cfc appearing in the Baez-Duarte formulation. It is shown that 

^ ■ R(x) can expressed by Ck and vice versa the sequence Ck can be obtained 

from the values of R(x) at integer arguments. We give also some rela- 
tions involving and R(x), in particular value of the alternating sum 
> ! of c fc . 



1. Introduction. 



The Riemann Hypothesis (RH) states that the nontrivial zeros of the function: 

1 ^ (-I)™" 1 



n=l 



where 3?(s) > and s ^ 1, are simply and have the real part equal to half, i.e. 
3?(s) = |. There are probably over 100 statements equivalent to RH, see eg. [T|, j2], 
In the beginning of XX century M. Riesz [3] has considered the function: 

R ^ = ]2 {k _ mm = x ^ kia2k + 2 y ( 2 ) 
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Unconditionally it can be proved that R(x) = 0(x 1/2+e ), see P §14.32. Riesz 
has proved that the Riemann Hypothesis is equivalent to slower increasing of the 
function R(x): 

RH R(x) = O (x l/4+£ ) . (3) 

A few years ago L. Baez-Duarte [Hj considered the sequence of numbers Ck defined 
by: 

He proved that RH is equivalent to the following rate of decreasing to zero of the 
above sequence: 

RH c k = 0(k~-* +t ) for each e > 0. (5) 

Furthermore, if e can be put zero, i.e. if Ck = 0{k~4), then the zeros of £(s) are 
simply. Baez-Duarte also proved in [S| that it is not possible to replace | by larger 
exponent. Although the title of the Baez-Duarte paper was A sequential Riesz-like 
criterion for the Riemann Hypothesis he has not pursued further relation between 
Ck and R(x). 

In this paper we are going to establish the relation between Ck and R(x). In 
Sect. 2 we will present formulae allowing to obtain values of R(x) and Ck much 
faster than from (j2j) and (jlj). In Sect. 3 we will use the fact that Ck can be obtained 
as forward differences of a appropriate sequence to express R(x) in terms of 
Next we will prove equivalence of the Riesz and Baez-Duarte criterion for RH. In 
the mathematical logic the iff obeys the transitivity rule: 

(p<=> q AND q <^> s) =>• (p <^> s) 

thus from © and © we have that R(x) = O (x 1/4+e ) c k = 0(k-'i +e ). 

However we will prove equivalence (Riesz criterion) (Baez-Duarte criterion) 
in a more general form, namely the exponents 1/4 and 3/4 will be replaced by 
arbitrary parameter 5 and combination 1 — 5: Ck = O (k~ 5 ) R(x) = O (x 1-5 ) . 
In the final Section we will speculate on some equations involving Ck and R(x), in 
particular we will calculate the alternating sum ^^L (— l) k Ck- 

2. Some facts on the R(x) and Ck 

The most comprehensive source of information about the Riesz function R(x) we 
have found on the Wikipedia For large negative x function R(x) tends to xe~ x . 
For positive x the behaviour of R(x) is much more difficult to reveal because the 
series (0) is very slowly convergent. Applying Kummer's acceleration convergence 
method gives 

fi{n) 
n 2 



R(x)=xJ2^re-^ (6) 



n=l 
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0.2 — i 




Fig.l The plot of R(x) for x £ (1, 20). Such a short interval is chosen to show the first 

zero of R(x). 
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Fig. 2 The plot of R(x) for x £ (0, 10 7 ). The part of R(x) smaller than -0.006 is skipped. 
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where \x is the Mobius function: 

1 if n = 1 

p(n) = ^ if n is divisible by a square of a prime (7) 
(— l) k if n is a product of k different primes 

Repeating Kummer's procedure gives: 

= +E^(«-*-0)- (8 > 

Using this formula we were able to produce the plot of for x up to 10 7 , see 
Fig.l and Fig.2. The first nontrivial zero of R(x) is xq = 1.1567116438.... The 
envelops on the Fig.2 (in red) are given by the equations 

y(x) = ±Ax*, (9) 

where A = 0.777506 . . . x 10~ 5 . 

It is very time consuming to calculate values of the sequence Ck directly from the 
definition (|1J), see [7], jH]. The point is, that for large j ((2j) is practically 1, and 
to distinguish it from 1 high precision calculations are needed. The experience of 
[H] showed that to calculate from (jlj) roughly Hog 10 (fc) digits accuracy is needed. 
However in jB] Baez-Duarte gave the explicit formula 1 for valid for large k: 



^ = Tk ^^gr + °w no) 

where the sum runs for nontrivial zeros p of C( s ) : CG°) = an d ^(p) 7^ 0. Maslanka 
in [7j gives the similar formula which contains the term hidden in o(l/k) in (jlOj) . 
Let us introduce the notation 

r(i - 2i) 

—— r -^— = a{p i )+ib{pi)=a i + ib i . (11) 

Assuming pi = \ + iji it can be shown jH] that and 6j very quickly decrease to 
zero: 



r i 



2 



C'(Pi) 



~ e~^ /4 (12) 



Finally we obtain for large fc: 

i ^ r / 7 aog(fc) ^ . ^tiogW M 

Cfe-i = -r 2^ S o< cos - bi sm V . (13) 



k ~ 4 i=i 



1 There is an error in |S] and there should be no minus sign in front of Ck—i in formulae (1.11), 
(1.12), (4.1), (4.11) in 0. 
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The above formula explains oscillations seen on the plots of published in jU] 
and 7], see Fig. 3 . Because these curves are perfect cosine-like graphs on the plots 
versus \og(k) it means that in fact in the above formula (|13|) it suffices to maintain 
only the first zero and skip all remaining terms in the sum. 



3. Relation between R(x) and cj. 

The values of Ck can be obtained as the first elements of the sequence of forward 
differences of the sequence: 

fO _ * f _ * f _ f _ fO _ \ 

Jo ~ /-(n\ ■'I — tlA\ Jl ~ tla\ J3 ~ /-/q\ J4 



C(2) Ji C(4) ' 2 C(6) Ji C(8) J4 C(10) 

(14) 

Then we form forward differences: 

/* = ft 1 ~ /w (15) 

and we have that = We will recall some facts from finite difference calculus 
adapted for our purposes [TUj: Let us define as usual the shift operator E: 



Ef(k) = f(k + 1). 
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Next we introduce sequence: 



Ck = (1 - E) fc /(0) = J2 ( ■) (-l)VO')- (16) 

3=0 



Then the following equalities holds [10J: 



e ^) /( 0) = ^|^, 

fc=0 



OO / \ JU 



e ,(i-^) /(0) = e - e -^ /(0 ) = e ^ /(*). 

fc=0 



from which it follows that: 

OO 



fc=0 fc=0 



In our case we put 

/(*) ' 



e(2fc+2) 

and finally we have: 

y c ± x k = e *y \zp = e _ R{x) . (17) 

^ k ^ k f(2k + 2) x K ' K ' 

Thus R{x) can be reconstructed from Ck- Vice versa, we will see later, see (}2T)j) . 
that within some accuracy Ck can be obtained from R(x). In the paper [S] it was 
suggested that the duality holds: 

Ck = (k- 5 ) <=> R{x) = O {x l ~ s ) . (18) 

Putting 5 = | — e gives original criteria Q and (JSJ). In fact we will prove it in the 
following form: 

Theorem 1. The sequence c k defined by decrease like Ck = O (k~ s ) if and only 
if the function R(x) defined by HJ) grows like R(x) = O (x 1 " 6 ) , wh ere 5 < 3/2. 

Remark: In fact 5 is smaller than 3/4, as shown by Baez-Duarte in jB]. 

Proof: The reasoning; that if c k = O (k- 5 ) then R(x) = O (x 1 ' 6 ) we will base 
on the following facts from Exercises 67 - 71 in the Part IV of famous book of G. 
Polya and G. Szego [TT]. We summarize these facts adapted for our purposes in the 
form: Let 

OO 

f(x) = Y,*kX k , (19) 

k=0 
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where are positive and decrease monotonically a > ai > a 2 . . . > > . . .. Let 
a be defined by 

i Hog A; 

loga fc (20) 

a 

and next the parameter b be determined from 

\ogf(x)~bx a . (21) 
Then for large x the following asymptotic relation is fulfilled: 

oo 

k5 akX k ~ (a6x Q ) 5 /(x). (22) 

fc=i 

In our case we have = l/k\ thus /(x) = e x and from Stirling formula we have 
a = 1 and next 6 = 1 and hence we have from above formula for large x: 

fc=i 

If we assume that IcJ < Afc -5 then we have 



^l£^ <AE -V (24) 



fc=l 



and from (|T7j) it follows: 

\R(x)\ < Ax 1 ' 5 (25) 

what is a desired inequality. 

We will show now the opposite implication: from R(x) = O (x 1-5 ) it follows 
that Cfc = O (k~ s )- In Appendix we prove the following inequality: 



R(k) 

Ck 



k 



< ^r^ 3/2 + ° (^ 2 ) (26) 

lo 



Because \ck\ — \R(k)/k\ < \ck — R(k)/k\ and we assume \R(k)\ < Bk 1 s thus we 
have 

\c k \ < Bk' 6 + O (VI) (27) 



To avoid nonsense 6 should be smaller than 3/2 and in fact Baez-Duarte showed jH] 
that existence of zeros on the critical axis requires 5 < 3/4. 

□ 



The comparison of the above bound (|26|) with real computer data is given in 
the Fig. 4. Here the fit (red line) was obtained by the least square method from 
the data with k > 10000 to avoid transient regime and it is given by the equation 
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y = 0.0117483x -1 ' 52655 . The fact that approximately Ck ~ R(k)/k was observed 
previously by S. Beltraminelli and D. Merlini [12] ■ It can be explained heuristically 
as follows: Baez-Duarte gives in [U] despite (J3J) a few formulae for c^. We need here 
the following expression being the transformation of (@J): 




For large k we can write: 




and comparing it with (JUJ) we get m R(k)/k for large k. 



10° -| 
10 _1 -| 
10" 2 - g 

10 _3 -B 



10^ -a 



10" 5 -i 
10" 6 - | 
10" 7 - | 
10" 8 -| 
10" 9 -i 
10' 



10" -i 
10 v 
10 ' :! 



10" 14 



m 

I 

3^ 



k - c * 



16 

least square fit 
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Fig. 4 The log-log plot of \c k - R{k)/k\ for k e (0, 10 6 ). 
4 Some other relations 

Using the formula (J2HJ) it is possible to calculate the alternating sum of c^: 



oo oo ^ ^ 



(30) 



A:=0 
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Numerically this sum is T l T=o(~ 1 ) kc k = 0.782527985325384234576688.... This 
number probably can not be expressed by other known constants, because the Simon 
Plouffe inverter failed to find any relation jT^j- By the Abel's summation the r.h.s. 
can be written as: 

1 1 r ( 1 \ C'(x) , , x 

Z»<m = 1+ J, (31) 

In fact more general than (|371|l formula holds: 

oo ^ oo / \ k , 



sj C(2fc + 2) 

where — 1 < s < ^ Here we have made use of the identity 



oo / 1 OO/ 

E s ' = T^^ (33) 

n=l v 7 fc=0 



The l.h.s. is convergent for — 1 < s < 1 while the r.h.s. converges for — oo < s < 1/2. 

The question of the convergence of the sum J2T=o Ck * s mucri more complicated. 
Formally summing both sides of ()16|) we get: 



J> = irVtf)) = /(-i) (34) 

fc=0 

As in our case f(k) = l/((2k + 2) we have 



because £(0) = — |, see e.g. 0, p. 19. The partial sums Y^l=o c k indeed initially 
tend from above to -2, but for n « 91000 the partial sum crosses -2 and around 
n ~ 100000 the partial sum starts to increase. These oscillations begins to repeat 
with growing amplitude around -2. The Fig. 5 shows the plot of distances of the 
partial sums J^ =0 from -2. Let us remark that at n ~ 10 8 the amplitude is rather 
very small: of the order 0.001. When we retain in (J 13)) only the first zero 71 it can 
be shown that this amplitude grows like n 1//4 , thus it appears that the above formal 
derivation (J53|) is wrong and the sum Ylk=o c k * s divergent. 
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Zc+2 



k=0 




20000000 40000000 60000000 80000000 100000000 

n 



Fig. 5 The distance from -2 of the partial sums Y^k=o f° r n = 1, . . . 10 8 . 

We have made analogous plot of partial sums Y^k=o(~ ^) k °k an d t nere we have 
seen oscillations around the limit value s = 0.782527985 ... of decreasing amplitude. 
Thus we speculate, that this partial sums behave as and R(x) accordingly: 



Et- 1 )' 



Cfc 



fc=0 



k=0 



On* 



(36) 



(37) 



Finally we would like to argue in favour of the two strange approximate equalities. 
Both follows from c& ~ R(k)/k for large k. The first follows when write this relation 
with the help of © and (@J): 



E 

i=o 



M(2j + 2) 



EC- 1 : 

3=0 



jj((2j + 2y 



(3f 



On both sides there appears inverses of C{2n). We have checked numerically that 
the difference between these two sums very quickly tends to zero. 

The second formula we get when in (fT7|) we put instead of Cfc simply R{k)/k: 



e 

—R{x) 

x 



J2m 



X 



k=0 



k k\ 



(39) 
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Thus we get R(x) as "entangled" combination of R(k) at positive integers. We end 
asking the question: Will such a kind of constraint help to prove (JHJ)? 

Acknowledgement We thank Prof. L. Baez-Duarte and Prof. K. Maslanka for 
e-mail exchange. To prepare data for some figures we have used the free package 
PARI/GP [Hj. 

Appendix 

In this appendix we will calculate the error of the approximation Ck ~ R(k)/k. 
Looking at (0) and (f2~KJ) we see that we have to estimate the sum: 



E 

n=l 



n 2 



e n 



k 



E 

n=l 



/i(n) 



n- 



n- 



oo 

<E 

n=l 



1 _ k 1 
— e 'r? = 



n- 



rr 



Instead of //(n) we have put 1. Let h(x) denote for 1 < x: 



(40) 



h{x) = — exp(-k/x ) 1 



ar 



x- 



This function is bounded by: 



n L/ , ^ 27 128 

< h(x) < — — ^TTT + 



2e 3 fc 2 ' e 4 k r 
and has one maximum. Thus we can apply the rule: 



°° PC 

h(n) < maxh(x) + / 

71=1 ^ 1 



h(x)dx. 



The integral is estimated as 



/ - (exp(-A;/a; 2 ) - (1 - l/x 2 ) k ) dx = / (e~^ - (1 - y 2 ) k )dy. 
Ji x Jo 



,- k y 2 d y = &r 1/2 



[\i-y 2 ) k dy 

Jo 



s/k 



4 fr 



-w 2 ^ < A;" 1 / 2 / e~ y2 dy 



2y/k' 



> 



ixk 



{ 2k )(2k + l) ~ 2k + 1 V 8A; 72k 



Here the Stirling formula in the form 

k\ = V2^k k k e- k+e{k \ 



12k + 1 

was used. Collecting all above estimations we obtain 

"OO J 



' < 9{k) < ' 



12Jfe 



;(exp(-A;/x 2 ) - (1 - l/x 2 ) k )dx < 



ar 



(41) 
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~ 2v^ 2A; + 1 V 8A; 72P / 16 144 
and finally from the starting sum (}4Ti|) we get the desired inequality: 

\R{k)/k - c k \ < ^#fc- 3/2 + %e^k~ 2 + ^V 5 / 2 + 128e~ 4 fc- 3 . 
16 2 144 

For fc > 16 it suffices to retain in the above inequality on the r.h.s only the leading 
term fc~ 3//2 . Let us remark that the integral ()41j) can be also taken from tables as it 
is the Euler Beta integral: 

and from 

B(a, x) ~ x~ a T(a) for x large 
see e.g. |15j §1 .8.7, we get: 

r(|) ^ 



f\l-y 2 ) k dy~ " l2 
Jo 



what for large k reproduces leading term in (}4"Tj) . 
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